Renormalization of the electron-spin-fluctuation interaction in the t — t' — U Hubbard 

model 
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We study the renormalization of the electron-spin-fluctuation (el-sp) vertex in a two-dimensional 
Hubbard model with nearest-neighbor (t) and next-nearest-neighbor {t') hopping by a Quantum- 
Monte-Carlo calculation. We distinguish between el-sp vertices involving interacting particles and 
quasiparticles, i.e. separate the renormalizations of the vertex from that of the quasiparticle residue 
1/Z. We show that for t' = 0, the renormalized el-sp vertex, not dressed by 1/Z, decreases with 
decreasing temperature at all momentum transfers. As a consequence, the effective pairing inter- 
action mediated by antiferromagnetic spin fluctuations is reduced due to vertex corrections. The 
inclusion of a flnite t' /t < 0, increases the Landau damping rate of spin fluctuations, especially in 
the overdoped region. The increased damping rate leads to smaller vertex corrections, in agreement 
with earlier diagrammatic calculations. Still, corrections reduce the spin-fermion vertex even at 
finite t' . 



I. INTRODUCTION 

The nature of the spin-fluctuation mediated interac- 
tion between charge carriers continues to attract high in- 
terest in the high-Tc superconductivity researchii^i^iii^i^. 
Theoretical calculations for the doped Hubbard and t— J 
models suggest that the strongest interaction between 
fermions is due to antiferromagnetic (AF) spin fluctu- 
ations. The low-energy description of this interaction 
has been advanced in the framework of the spin-fermion 
model^. A magnetically mediated interaction has been 
intensively studied since the early days of high-Tc era 
because AF spin fluctuations induce d— wave pairing be- 
tween fermions Various numerical and analytical cal- 
culations based on Hubbard and t — J models, as well as 
semi-phenomenological spin-fermion models, have been 
performed in the normal and superconducting states of 
the cuprates. The results of these calculations are in 
quite good agreement with a large number of experimen- 
tal data. In particular, these calculations show that the 
magnetic response in a d— wave superconducting state 
contains a resonance peak, in addition to a gapped con- 
tinuum. The location of the peak, and its "negative" 
dispersion are in agreement with experimental data. It 
has been argued- that the scattering from the resonance 
is strong enough to explain the pcak-dip-hump feature 
in the fermionic spectral function, the 5— shape disper- 
sion for antinodal fermions, and similar features in SIS 
tunneling, optical conductivity, and Raman response. 

The strength of the feedback effects from the reso- 
nance peak on fermions, and the magnitude of the spin- 
mediated Tc depend on the size of the spin-fermion cou- 
pling. The estimates for the coupling strength still vary 
substantially. Abanov et al argued^ that the coupling 
g is comparable to the Hubbard U and is quite strong, 
g ^ 0.7eV. On the other hand, Kee et al^ argued 
that the effective spin-fermion coupling is much weaker 



g ~ O.OlieV. In the latter worl^, the coupling was ex- 
tracted from the specific heat data. To a large extent, 
the difference between the two results is related to dif- 
ferent choices for the fermionic density of states Nq: a 
large value of g is obtained by assuming a large Lut- 
tinger Fermi surface, with Nq ~ leV~^ (Ref. '?'), while 
a smaller g is obtained by assuming that the density of 
states Nq ^ ^ 10eV~^ is the same as in a weakly 
doped quantum AF. 

This difference brings in the issue of how much the full 
spin-fermion vertex gtq differs from the bare gj^^ due to 
vertex corrections. Here, k and k + q are incoming and 
outgoing fermionic momenta, and q is the bosonic mo- 
mentum. The ratio r(k, q) = .9kq/3kq determines the 
vertex renormalization originating from electronic corre- 
lations. Some of us have recently shown numerically that 
this quantity is substantially renormalized by strong elec- 
tronic correlations both for the electron-phonon (el-ph) 
verte3*i^ as well as for the el-sp vertex™. 

Within the spin-fermion model, the role of the vertex 
renormalization is well understood in the limits of very 
small and relatively large doping. At very small dop- 
ing X, when long-range AF order is still present, nearly 
all carriers are localized, and doped fermions form small 
pockets around (7r/2, 7r/2) and symmetry related points. 
In this situation, the full spin-fermion vertex g^q vanishes 
when q coincides with the antiferromagnetic momentum 
Q = (7r,7r). The vanishing of g^Q is the consequence of 
the Adler principle that the interaction should preserve 
massless Goldstone bosons in the theory and, therefore, 
should vanish at the ordering momentum of Goldstone 
bosons. For q 7^ Q, the vertex is non-zero, but small in 
|q — Q|/|Q|. Mathematically, the vanishing of g^Q is the 
result of dressing up of the bare interaction g^q by coher- 
ence factors, associated with the antiferromagnetic order. 
From this perspective, the strong reduction of g^q when 
q is close to Q is the result of a strong vertex renor- 
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malization. SchriefFer arguedi^, that, if the pocket-hke 
Fermi surface survives even when long-range AF order is 
lost, the strong vertex renormalization extends into the 
paramagnetic phase. According to Ref. [H 



r2(k,q)cx [(q-Q)^ 
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[X(q,^ = 0)]' 



(1) 



where ^ is the AF correlation length. Once Eq. ([T]) is 
valid, the effective spin-mediated interaction 



F(k,q)cx|r(k,q)|2x(q), 



(2) 



is considerably weaker than without vertex renormal- 
ization, and, most importantly, it is no longer peaked 
at q ^ Q. This is an important fact, since the peak 
structure at q ~ Q was responsible for a strong attrac- 
tive d— wave component in the pairing potential. The 
magnetically-mediated d— wave pairing is still possible 
even in this situation, but the corresponding Tc is much 
smaller than one would get without vertex correctionsi^. 

In the opposite limit of large doping, fermions are 
itinerant, and have a large Fermi surface which crosses 
the magnetic Brillouin zone boundary at hot spots (this 
requires a sufficiently strong ncxt-ncarest-neighbor hop- 
ping). A spin- fluctuation with momentum near Q then 
can decay into a particle and a hole. Analytical calcu- 
lations within the spin-fermion model show that in this 
situation, the vertex is only weakly renormalized, and 
actually increases with respect to its bare value, i. e. 
5kQ > 9kcr^- Moi^s precisely, the relation between (7kQ 
and g'^Q depends on the angle (j) < tt between Fermi ve- 
locities at hot spots separated by Q^: 



9kQ 



1 



(3) 



where is of order 1 . Equation ([3|) shows that the renor- 
malization is largest when (f) n. This limit corresponds 
to an almost nested Fermi surface at hot spots. In op- 
timally doped cuprates, the velocities at the hot spots 
separated by Q are almost perpendicular to each other, 
i. e., (p ~ 7r/2. For ^ ~ 2, the vertex correction is then 
about 10% (even smaller, ^ 4%, when one adds regular 
termsi^J^). 

The difference between the itinerant and localized lim- 
its raises the question how vertex corrections evolve as 
the system approaches half-filling. One attempt to ad- 
dress this issue was undertaken in Ref. [l^. The authors of 
Ref. 16 considered a toy spin-fermion model in which one 
can study the evolution from the Luttinger Fermi surface 
towards hole pockets. This evolution involves a topolog- 
ical transition below which a pocket is splitted from a 
large Fermi surface, and the Luttinger theorem breaks 
down. Within this model, vertex corrections evolve to- 
gether with the Fermi surface: they are small and positive 
(i- e. (7kQ > Skq) i^i the limit when the Fermi surface is 
large, but change their sign near the topological transi- 
tion, become negative and rapidly increase in magnitude 



as the Fermi surface evolves towards hole pockets. For a 
pocket-like Fermi surface, vertex corrections almost can- 
cel the bare vertex, and the full gkq coincides with Schri- 
effer's result, Eq. ((J). 

In this paper, we address the issue of the strength 
of vertex corrections in the two-dimensional t — t' — U 
Hubbard model by means of a quite accurate numerical 
method - the determinantal Monte-Carlo^- algorithm. 
The purpose of our calculation is to explore the extent of 
validity in doping and temperature of Schrieffer's argu- 
ment^i^ as well as of the spin-fermion model calculation^: 
Whereas the analytical results are only valid when the 
correlation length ^ is large, Monte-Carlo calculations 
are not restricted by this condition. We also want to get 
an insight as up to what extent a description in terms 
of low-energy spin degrees of freedom, i. e. in terms 
of a spin-fermion model is appropriate for the Hubbard 
model. 

Vertex corrections have been analyzed by some of us 
in Ref. [HI The present paper extends this earlier letter 
in three different directions. First, in order to separate 
vertex and quasiparticle renormalization effects, we ana- 
lyze the vertex renormalization and the pairing interac- 
tion with and without the inclusion of the wave-function 
renormalization Z . The vertex without the renormaliza- 
tion of Z is labeled by r(p, q) before and in Eq. ([7]) below. 
This vertex is for the interaction between spin fluctua- 
tions and strongly interacting particles. This vertex has 
to be used for comparison with analytical results which 
don't include the rescaling by the fermionic Z. The fully 
renormalized vertex is the one rescaled by Z (we label it 
by 'j{p,q), see Eq. Q below}!^. It describes the inter- 
action between spin fluctuations and quasiparticles. We 
discuss the difference between the two vertices in Sec. lIIIi 
and also discuss there the behavior of the pairing inter- 
action between interacting particles and quasiparticles. 
Second, we calculate the momentum dependence of the 
vertex. Third, we discuss in detail the effect of bosonic 
damping on the strength of the vertex renormalization. 
We show that the introduction of a next-nearest-neighbor 
hopping t', which makes the bosonic damping possible, 
substantially reduces the renormalization of F (but not 
of 7). 

Our numerical Monte Carlo simulations are performed 
on a 8 X 8 lattice at different doping densities and different 
temperatures. In our flnite lattice, the (tt, 0) point is the 
one closest to the hot spots, so that charge carriers near 
the (tt, 0) region are strongly affected by antiferromag- 
netic spin fluctuations. Therefore, we will mainly con- 
centrate on the particular scattering process in which the 
incoming electron and spin fluctuation carry momenta 
p — (— TT, 0) and q — (tt, tt), respectively. Within our p- 
points mesh, the points p and p -I- q lie sufficiently close 
to the Fermi surface. 

Our paper is organized as follows: In Section |TT1 we 
define the Hamiltonian and describe the numerical ap- 
proach, which is based on the Quantum-Monte-Carlo 
evaluation of the linear response to an external spin per- 



3 



turbation. In Section IIIIl we present our numerical re- 
sults and discuss their qualitative relation with the spin- 
fermion-model calculation and Eq. [1] Finally, in Sec- 
tion |IVl we discuss in detail our main conclusions. 



II. 



MODEL AND NUMERICAL APPROACH 



Our starting point is the one-band Hubbard model, 
H = ^ (cL V + 4q J - i' ^ {clc^a + c]^c,,) 

+U'^n,inii, (4) 

i 

Here, the operators cf^. and c^^ as usual create and de- 
stroy an electron with spin a at site i, respectively, and 
(ij) and {{ij)) denote a sum over nearest and next- 
nearest neighbor lattice sites i and j, respectively, riia- — 
^\a^ia number operator. Finally, U is the onsite 

Coulomb interaction and the nearest-neighbor hopping t 
is chosen as the unit of energy. In order to consider the 
effects of damping, we include a next-nearest-neighbor 
hopping term t' . 

In our simulation, we use a linear-response approach 
(see also Ref. 10) in order to extract the el-sp vertex. 
In this method, one formally adds to Eq. (|4|) the interac- 
tion with a momentum- and (imaginary) time-dependent 
spin-fluctuation field in the z-direction S'qC"**"'^— in the 
forni-^^ 



H, 



el— sp 



kq(T 



(5) 



where is the bare el-sp coupling (equal to the Hub- 
bard U in the one-band Hubbard model) . One then con- 
siders the "anomalous" single-particle propagator in the 
presence of this perturbation defined asi^ 

G^(p,g) EE - rdre^(P"+'"')-(T.Cp+q.(T)c^,(0))ff+ff„,., 
Jo 

(6) 

Here Qh+h^i^^-p is the Green's function evaluated with 
the Hamiltonian H + Hci-sp- Diagrammatically, GAip, q) 
has the structure shown in Fig.[l] so that the el-sp vertex 
r(p, q) can be expressed quite generally in terms of Ga 
and of the single-particle Green's function G{p) in the 
form 

1 1 1 GAip^q) 



r(p, <?) = lim, 



°9LS^1 + Uxzz{q) G{p + q)G{p) ' 



(7) 



with Xzz (q) the longitudinal spin susceptibility. Because 
only the limit — > is relevant in Eq. [7l it is sufficient 
to calculate the leading linear response of Ga to iJoi-sp, 
which is given by 

r/3 



GAip,q) = S^ I dTe^(Po+,o)r / dr'e-^«"^'^.g;jqX 
Jo Jo 





G.(p,q)= 



FIG. 1: Diagrammatic representation of Ga{j>,<i) within lin- 
ear response to Sq. The thick solid lines represent dressed 
single-particle Green's functions of the Hubbard model. The 
wavy line denotes the external perturbation in Eq. ([5}. The 
dashed line represent the Hubbard interaction U and the black 
ellipse stands for the longitudinal spin susceptibility X^^il)- 



where 0"*" is a positive infinitesimal. Notice that Sq can- 
cels in Eq. [T] The two-particle Green's function in Eq. dH) 
is evaluated with respect to the pure Hubbard Hamilto- 
nian (Eq. (g])). 

Within a quasiparticle approach, one can shift the 
wave-function renormalization factor Z{p) > 1 from 
the Green's function into the definition of the effective 
el-sp vertexfiS Thus, for states close to the Fermi en- 
ergy with a single-particle Green's function of the form 
G{p) = Z{p)^^ {i pq — Ep)~^, where Ep is the quasi- 
particle energy, the effective quasiparticle-spin coupling 
is given by: 



lip, q) 



np,q) 



^Z{p)Z{p + q) 



(9) 



Z{P) 
vertex 



be- 



Numerically, Z{p) is evaluated as 
Im[l/G(p)]/p(^ii^. Therefore, 7 is the 
tween quasiparticles and spin fluctuations. 

The total pairing interaction Vp for the exchange of a 
single spin fluctuation can be expressed in terms of the 
vertex F as 



Vp{p,q)^\Tip,q)\^ -U^ -Xzziq) 
where Xzz{q) is the spin susceptibility 
1 



(10) 



Xzziq) 



dr 



and 



(T.m^(r)m^_q(0)) , 

^ II . (11) 



'J 



Eq. [TO] describes the effective pairing interaction between 
interacting particles. Including the wave-function renor- 
malization as in Eq. [51 we introduce the effective pair- 
ing interaction between quasiparticles (this is in complete 
analogy to the treatment of the strong-coupling super- 
conductors in Refi^) 

Vp{p, q) 



,(r' + 0+)c,.,(r')cp+q.(r)ct,(0))H,(8) ^^^^^^ = Z{p) z[p ^ q) = l^^^''^)' ^ ^^^^^ ' ^^'^ 
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(0,0) {7t/2,n/2) (7i,7t) 
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FIG. 2: (color online) (a) Spin susceptibility Xzzil)^ (b) the 
renormalized el-sp vertex |r(p, g)|, and (c) the total pairing 
interaction Vp{p,q) as a function of spin-fluctuation momen- 
tum transfer q. Here U = 4, t' — 0, p = (— 7r,0), and the 
doping density S = 0.12. The value of the inverse temperature 
/3 is indicated by the shape of the symbol. 



Clearly, both 9) and Vp{p,q) contain the contri- 

butions from both the real and imaginary parts of the 
vertex. Within a Quantum-Monte-Carlo approach, the 
effects and strength of the particle-particle vertex in a 
given channel can be also indirectly measured by compar- 
ing the full pairing susceptibility with the corresponding 
"bubble" approximationSi. 



III. RESULTS AND DISCUSSIONS 

In order to explore the momentum structure of the el- 
sp interaction, we first plot Xzziq), |r(p, q)\, and Vp{p, q) 
versus q at different temperatures in Fig. [21 Here, the 
spin- fluctuation momentum transfer q is along the (1, 1) 
direction and U is at an intermediate coupling value, i.e. 
U ~ A. One can readily see that both Xzz{q) and Vp{p, q) 
are peaked at momentum transfers around the antifer- 
romagnetic vector Q = (tt, tt), and the strength of the 
peaks increases when the temperature is lowered. This 
demonstrates that the q- and T-dependences of the spin 
susceptibility x dominate the temperature behavior of 
the pairing interaction, despite the reduction of the ver- 
tex r as the temperature is decreased. From Fig. [2l we 
observe that the decrease with temperature of T at large 
momentum transfers is stronger than at small momentum 
transfers, indicating that the vertex correction at q ~ Q 
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FIG. 3: (color online) The absolute values of the (complex) 
renormalized el-sp vertices |r(p, q) \ and |7(p, 5)! vs T at [/ = 4 
and ?7 = 8 at the doping density 5 — 0.12. We set p — 
(— 7r,0), q — (tTjIv), and t' = 0. 



is larger than at small q. This q-dependence of the el-sp 
vertex is qualitatively in agreement with the prediction 
of Eq.([T]). Our finding is also in good agreement with 
the work of Bulut et al^ (carried out at U = 4), which 
shows a value of g = 0.8 (corresponding |r| = 0.8 in our 
notation) can produce an effective coupling glJ which is 
consistent with the results of Monte Carlo calculation 
of the irreducible particle-hole vertex, and that the ef- 
fective particle-particle interaction originating from the 
Hubbard U increases with lowering temperature and can 
reach large values. 

Monte Carlo results for T{p,q), '){p,q), Z{p), Vp{p,q) 
and Vp{p,q) are displayed in Figs. [3] and S] for intermedi- 
ate as well as strong correlation. We compute the vertices 
at zero bosonic Matsubara frequency and the smallest 
fermionic Matsubara frequency uji = ttT. Because loi is 
finite, and also because at a nonzero doping, the exci- 
tation spectrum is particle-hole asymmetric, the vertices 
r and 7 have both real and imaginary parts (the imagi- 
nary part obviously vanishes at T = 0). We notice that 
in the underdoped regime {6 — 0.12) for both interme- 
diate {U — 4) and strong correlation {U — 8), both |r| 
and I7I are strongly renormalized below a characteristic 
temperature {T < J = 0.5 for U = 8). Although our 
simulation cannot reach low temperatures, a clear trend 
is observable in both |r| and I7I, which tend to go to 
small values at low temperatures, at least for U = 8. A 
more careful look shows that the reduction of the ver- 
tices r and 7 is chiefly due to the reduction of ReT and 
Rej. The imaginary parts of the vertices are small at 
weak and intermediate correlations (at our lowest acces- 
sible T), but for ?7 = 8 they can become comparable to 
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FIG. 4: (color online) The pairing interaction Vp{p,q) vs T 
for (a) (7 = 4 and (h) U = 8 at the doping density 5 = 0.12 
(full symbols, the solid line is a guide to the eye). The open 
symbols (dashed line) show the RPA results. In panel (c), we 
show the T-dependence of the wave-function renormalization 
Z{p) and in panel (d) the effective pairing interactions Vp{p, q) 
between quasiparticles. In all panels, p — {—it, 0), q = (vr, tt), 
and t' = 0. 



the real parts. 

Further, our numerical results presented in Fig. 3] 
clearly show that, in both intermediate- and strong- 
correlation regimes, the renormalized pairing interaction 
Vp (which does not include the dressing by the quasipar- 
ticle Z) is smaller than the RPA result obtained with the 
full susceptibility Xzz- The difference becomes strongest 
at the lowest accessible temperatures. The qualitative 
temperature behavior of the effective pairing interaction 
between quasiparticles Up, is quite similar to Vp in the 
intermediate-correlation {U = 4) regime. However, it is 
quite different in the strong-correlation ([/ = 8) regime, 
where it displays a mild increase or a saturation at low 
T's. Based on the results shown in Fig. 21 we conclude 
that for t — 0, both vertex corrections and the renormal- 
ization of Z are important for the spin-mediated pairing 
interaction. Vertex corrections suppress Vp (p, q) com- 
pared to the case of free fermions. The rescaling by Z fur- 
ther suppresses the effective pairing interaction between 
quasiparticles, Vp{p,q). Note, though, that the pairing 
interaction Vp{p, q) at a bosonic q — (tt, tt) and fermionic 
p near the antinodal point still increases with decreasing 
temperature. 

Fig. 13] shows that vertex corrections are substantial 
even in the case when the AF correlation length ^ is quite 
small^?. i. e. of order of the Cu-Cu distance. This is in 
contrast to the situation discussed by Schrieffer in which 
holes move in an AF background, which is unaffected by 



3 
cr 



1 1 1 
— • — U=4 ^^^^^^^ 


1 1 


— ^U=8 




1 1 1 







0.0 0.2 0.4 0.6 O.i 
T 



1.0 



FIG. 5: (color online) Vj,{p, q) /Z{p) vs T for [/ = 4 and t/ = 8 
at the doping density 5 = 0.12 (full symbols, the solid line is 
a guide to the eye), p and q are the same as in Fig. |4] 



the charge carriers, i. e. ^ is large. How can we under- 
stand this result? In the AF precursor (large-^) case^^^, 
the vanishing of the el-sp interaction in the long-range 
ordered AF state is the result of dressing up of the bare 
interaction by the AF coherence factor, which is small 
at the top of the valence band. The coherence factor is 
due to the interference effect of the quasiparticle, which 
forms a "spin-bag" , i. e. a hole dressed by a short-range 
AF background. Although in our calculation the AF pre- 
cursor is no longer present, it is well established from our 
earlier QMC work on the evolution of the single-particle 
spectral function A(k, from an insulator to a metal 
that below T ^ J, the electronic excitations display an 
essentially doping-independent feature. More precisely, a 
"band" of width J forms, in which "spin-bag" -like quasi- 
particles propagate coherently. The continuous evolu- 
tion can be traced back to one and the same many-body 
origin: the doping-dependent AF spin-spin correlation. 
Therefore, we argue that the interference effect of the 
spin-bag quasiparticle plays a similar role in reducing 
the el-sp vertex in the strongly-correlated undcrdoped 
regime. 

Above we considered the pairing vertices for interact- 
ing particles and for quasiparticles. The actual strength 
of the pairing interaction is proportional to the product of 
the effective interaction Vp and the density of states at the 
Fermi level N{Ep). Since we restricted with only the pole 
component of G, the conservation of the particle num- 
ber implies that N{Ef) is increased by a factor Z with 
respect to the noninteracting density of states No{Ep), 
i.e., N{Ef)vp = No{Ef) Vp/Z. This is similar to the 
case of the electron-phonon mediated pairing, where the 
interaction is also rescaled by one power of Z (this last 
rescaling is a well-known McMillan result A A/(l -I- A). 
Our consideration, in which we defined Z as the overall 
wave-function renormalization is formally different from 
electron-phonon problem, where the fermionic self-energy 
depends only on frequency. However, the final result is 
the same as in the case where self-energy only depends 
on frequency, the density of states is not renormalized 
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FIG. 6: (color online) (a) Spin susceptibility x^^il)! (b) 
"damping rate" jd, and (c) el-sp vertex |r(p, q)\ as a function 
of doping density S for U — 4. The value of t' is indicated by 
the shape of the symbol. Here p = (— 7r,0), q = (vr,7r), and 
/3 = 4. 



by Z, but one power of Z, compared to Vp ^ Vp/Z'^ is 
eliminated by the mass renormalization m* /m = Z. 

In order to have a measure of the pairing strength 
including the renormaUzation of the density of states, 
we thus plot the quantity Vp{p,q)/Z{p) in Fig. [5] as a 
function of temperature for 6 = 0.12. We observe that 
Vp{p, q)/Z{p) increases with decreasing temperature in 
both intermediate- and strong-correlation regimes. On 
the other hand, the suppression of Vp due to vertex cor- 
rections and the increase of Z with decreasing tempera- 
ture considerably reduce the pairing (i. e., without vertex 
corrections) , in particular in the low-temperature regime. 

As discussed in the Introduction, previous work on 
the spin-fermion model^ suggests that the vertex correc- 
tion |r — 1| (not including Z) gets considerably reduced, 
whenever spin fluctuations get damped. To explore this 
fact, we have compared numerical results for different 
values of the next-nearest-neighbor hopping {t' — and 
t' = —0.25), since one expects the damping to increase 
for larger negative t' . The spin susceptibility Xzz, a quan- 
tity proportional to the damping rate 7^, and the el-sp 
vertex F are plotted in Figs. [Slja)-[nijc) as a function of 
doping density. When t' — 0, the Fermi surface is par- 
ticlelike and encloses the zone center F = (0, 0). On the 
other hand, when t' = —0.25, the Fermi surface is hole- 
like and encloses the zone corner Q = (it, tt), and, in par- 
ticular, it crosses the magnetic Brillouin Zone boundary 
at so-called hot spots. In this second case, spin fluctu- 



ations get strongly damped due to the decay into parti- 
cle and holes near the Fermi surface. Due to the finite 
temperature and short correlation length, a damping of 
spin fluctuations is also present for t' — 0. The quantity 
7d = b/c shown in Fig. [6l is obtained by fitting the spin 
susceptibility x(g,iwm)~^ to the form a + bujm + cw^, 
and is, thus, proportional to the spin damping rate. In 
the underdoped region, the spin susceptibility is strongly 
suppressed by a negative t' due to its frustrating effect on 
the AF alignment, whereas the difference of jd is small, 
as seen in Fig. [D^a) and Fig. ^h). On the other hand, 
in the overdoped region, where the holelike Fermi sur- 
face passes close to (tt, 0) and symmetry points (whereas 
the particlelike Fermi surface moves away from the mag- 
netic Brillouin zone boundary), t' has little effect on the 
spin susceptibility, while the spin damping rate is dra- 
matically increased by a negative t'. Fig. ^c) clearly 
shows that the el-sp vertex |r| is larger for t' = —0.25 
than for t' = 0, demonstrating that the vertex correction 
is reduced either by suppressing the spin susceptibility 
or by increasing the spin damping rate. However, the 
largest difference in |r| is not observed at the smallest 
doping density S = 0.066, where the suppression of the 
spin susceptibility is strongest, or near the doping den- 
sity 5 — 0.20, where the difference of the damping rates 
is largest. This result suggests that, although both the 
spin susceptibility and the damping of spin fluctuations 
are related to the magnitude of vertex corrections, they 
are probably not the only relevant factors. 

In Figs.[7I^a) and [7Kb) we show the effects of the next- 
nearest-neighbor hopping t' on the effective el-sp vertex 
7(p, q) and the effective pairing interaction Vp{p, q)/Z{p). 
In contrast to F, the vertex 7 ~ T / Z shows a weak depen- 
dence on t' (and also on 5), and in the overdoped regime 
is even smaller for t' — —0.25 than for t' = 0. The differ- 
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ent behavior of the vertices T and 7 is obviously due to 
the Z factor, which increases with increasing \ f\ and/or 
decreasing doping. The result fo Vp{p, q)/Z{p) shows that 
the pairing strength has almost no dependence on t' in 
the optimal and ovcrdoped regimes, and is somewhat re- 
duced at a negative t' in the underdoped regime. 

Finally, we comment on the overall sign of the vertex 
correction. In our cases studied, as shown in Fig. [2] and 
Fig- El |r| is always less than 1, indicating that the el-sp 
interaction is suppressed by vertex corrections. This is in 
contrast to the finding in the spin-fermion model, which 
shows that the total vertex correction at moderate dop- 
ing has a positive sign, i. e., vertex correction actually 
increases the el-sp interaction. This effect, however, is 
rather small and for moderate ^ may be overshadowed 
by the contributions from high-energies {0{Ef)), not in- 
cluded in the low-energy spin-fermion models. The com- 
parison between our calculations and the results in the 
spin-fermion model shows that, while the effect of damp- 
ing of spin fluctuations on vertex corrections is quite ro- 
bust and agrees with low-energy considerations, the sign 
and the magnitude of the vertex correction near optimal 
doping may be model dependent. 

IV. CONCLUSIONS 

In summary, based on Quantum-Monte-Carlo simula- 
tions, we have studied the renormalization of the el-sp 
vertex in the two-dimensional t — t'~U Hubbard model. 
We found that the fully renormalized el-sp vertex be- 
tween bosons and free fermions decreases quite gener- 
ally with decreasing temperature at all spin-fluctuation 
momentum transfers. We distinguished between vertex 
renormalization itself and the renormalization due to the 



dressing of the vertex by quasiparticle Z factors. We 
analyzed the effect of adding a negative next-nearest- 
neighbor hopping term t' to the dispersion. This term 
changes the topology of the Fermi surface and allows 
a spin fluctuations to decay into an electron- hole pair. 
We found that t' reduces vertex corrections (but not Z), 
in agreement with previous results on the spin-fermion 
model. However, in contrast to Ref. we did not ob- 
serve a positive vertex correction, i. e. the renormalized 
vertex is always smaller than the bare one. 

The suppression of the fully renormalized el-sp ver- 
tex, particularly at t' = 0, gives rise to a substantial 
reduction of the pairing mediated by antiferromagnetic 
spin fluctuations in both the intermediate- and strong- 
correlation regimes. This result extends Schrieffer's ar- 
gument (Eq. [T|) about the suppression of the el-sp vertex 
to the case in which the AF precursor is no longer present. 
Notice, however, that F and x behave in the opposite way 
(as implied by Eq. [l] only when one considers their tem- 
perature behavior. In contrast, when doping decreases at 
a fixed temperature, F increases with increasing x (the 
results are not shown here). Thus, in the situation when 
the system is away from the AF precursor, there is no 
general proportionality relation implied in Eq. [TJ 
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